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CILAPiER I 



Summary 

Recovery factor and friction factor data are presented here for air at 
a stagnation temperature of 1 1 0°F flowing through a cylindrical tube at 
Mach numbers ranging from 1.3 to 2.6 and at inlet diameter Reynolds 
numbers from 0.47 x 10^ to 4.8 x 10^. 

Varying from a value of about 0.81 at a Mach number of 1.3 to about 
0.905 at a Mach number of 2.6, the recovery factor has been shown to be 
an almost linear function of only Mach number over the entire range of 
inlet diameter Reynolds numbers. Although the results obtained indicate 
three flow regimes, there appears to be no dependency of the recovery 
factor on the diameter Reynolds number at inlet. In addition, there 
apparently is no simple relationship between recovery factors and Prandtl 
number involving fractional exponents. 

The mean apparent friction factors are in good agreement with the 
results of Keenan and Neumann^ ^ and Ketchum^^. From a plot of mean 
apparent friction factor as a function of length Reynolds number three 
flow regimes which are essentially laminar, transitional, and turbulent 
can be identified with the transitional regime apparently lying between 
length Reynolds numbers of 15 x 10^ and 30 x 10^. The friction factor 
attains a minimum value of about 0.0013 at the transition region and 
after transition decreases steadily with increasing length Reynolds 



number from a maximum value of about 0.0035. 
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Preliminary results were obtained with the recovery factor apparatus 
modified to reduce heat transfer to the test section;- although neither ex- 
tensive nor conclusive, these data indicate that the recovery factors 
obtained will be from one-half to one percent lower than those determined 
with the unmodified apparatus and that recovery factor is independent of 



stagnation temperature. 



CHAPTER II 



Nomenclature 

Symbol Meaning 

A cross-sectional area of test section (in^) 

A* cross-sectional area of nozzle throat (in-) 

Cp specific heat at constant pressure (Btu/lb-F) 

Cw nozzle flow coefficient (dimensionless) 

D diameter of test section (in) 

D'-^ diameter of nozzle throat (in) 

f local apparent friction factor; friction factor which is the mean 

of the values for two adjacent stations in the test section 
(dimensionless) 

f mean apparent friction factor; friction factor which is the 

mean of the values for the first station and another station in 
the test section (dimensionless) 
g gravitational acceleration (ft/ sec^) 

G mass velocity in test section (lb/ sec-ft^) 

mass velocity in nozzle throat (ib/sec-ft^) 
hj height of mercury in leg of McLeod gage (cm) 

h2 height of mercury in capillary tube of McLeod gage (cm) 

J proportionality constant (ft-lb/Btu) 

k specific heat ratio (dimensionless) 

k thermal conductivity (Btu/hr-ft-F) 



Symbol 



Kleanin^q 



^^max 

M 

p. p 



Pr 

R 

R’ 

Re 



Re 



- s 



longitudinal distance from exit plane of nozzle to a point in the 
test section (in or ft) 

longitudinal distance from entrance to cylindrical tube (exit 
plane of nozzle) to a point where the Mach number is 1.0 (in or ft) 
Mach number (dimensionless) 

static pressure; P^ar “ barometer; Fjs " downstream stagnation 
tank; Pgas- gasometer tank; Pd- upstream of orifice;AP(j- 
pressure differential across orifice; upstream stagnation 

tank; P^,- wall of test section; Pvac” evacuated chamber 
surrounding test section (psi or psia) 

Prandtl number (dimensionless) 
gas constant (ft -lb/ lb -F) 
recovery factor (dimensionless) 

Reynolds number; Rej^j- based on test section diameter; 

Rej^- based on distance from exit plane of nozzle; Re^- based 
on distance from nozzle Jthroat (dimensionless); Rej^*!'- based 
on nozzle throat diameter 

mean Reynolds number; value v/hich is a mean for two adjacent 

stations of test section (dimensionless) 

radius of station temperature thermocouple (in) 

radial distance between differential temperature thermocouples 



Symbol Meaning 



t 



T 

w 



X 



0 



radius of test section wall (in) 

temperature (^F); t^w- adiabatic wall temperature; down- 
stream stagnation temperature; temperature in gasometer; 

tj^- mean stream temperature; room temperature; 

t - temperature at station thermocouple; ^ t^ - temperature 
differential between differential temperature thermocouples; 
t^g- upstream stagnation temperature; 
temperature (°Fabs); same subscripts as for t 
mass rate of flow (ib/sec); Wg- mass rate of flow through 
orifice; Ws- isentropic mass rate of flow; w=!=- mass rate of 
flow .through nozzle; w^^g- mass rate of flow into gasometer 
longitudinal distance from nozzle throat to a point in the test 
section (in) 

viscosity (ib/sec-ft); yUm - viscosity evaluated at mean stream 

temperature; viscosity at nozzle throat 

4fL 

dimensionless group, rnax 
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CH^FTER III 



Introduction 

The experimental work reported in this thesis is part of Project 
Die 6418 sponsored by the Office of Naval Research to determine recovery 
factors, friction factors and heat transfer coefficients for air flowing at 
supersonic velocities, NIore extensive kiiowledge of the heat transfer 
coefficient is becoming increasingly important in the design of aircraft, 
aircraft coolers, missiles, compressor and turbine blading. However, the 
heat transfer coefficients as defined for air flowing at supersonic speeds 
cannot be determined without a prior knowledge of recovery factors, data 
for which is currently insufficient and inadequate. It is the purpose of 
this thesis to determine the magnitude of recovery factors and the 
parameters upon which the recovery factor depend in the flow of air at 
supersonic speeds through a cylindrical tube. In addition, data on closely 
related friction factors are to be obtained incidental to the determination 
of recovery factors. 



Recovery Factor 

The coefficient of heat transfer for compressible flow is defined as 
the rate of heat flow per unit area per degree temperature difference 
between wall temperature and the adiabatic v/all temperature: 



h = 

The temperature difference » Tw 



q/A 
T\v - T 

- Taw* is 



(1 

aw 

now almost universally accepted 



where any compressibility effects are encountered since its use yields 



values of h which are independent of the temperature difference at high 
subsonic and supersonic speeds. In earlier investigations, however, the 
heat transfer coefficients were often based on the temperature difference 
between and Tm or Tq. 

Sine 0 ^aw is defined as the temperature the wall would have if it 
were in equilibrium with the stream and thermally isolated from its 
surroundings, i.e., an adiabatic wail, it carmot be measured when heat is 
being transferred by the wall. However, by determining the value of the 
recovery factor defined in terms of Taw &nd the stream properties found in 
an adiabatic test. Taw can be predicted under the conditions when heat is 
being transferred by the wall. 

The definition of the recovery factor used is the dimensionless 
quantity; 



R* 



T - T 
^ aw ^ m 



T - T 



( 2 ) 



in which the stagnation temperature can be measured at a point in the 

stream where the velocity is zero; can be measured in an adiabatic 

test; and the mean stream temperature can be evaluated from a one- 
dimensional analysis of v/all and stagnation pressure measurements. 

Since in steady, one dimensional flow of a perfect gas in a duct of 
constant cross-sectional area in the absence of heat exchange 



T - T 
^ o m 






( 3 ) 



where V is the velocity of the mean stream at the point where Tj^ is 
evaluated, 



R’ = 



aw 



- T 



m 



vV 



(*) 



and 



Taw = T^ + R* 



V2 



(5) 



Equation (5) thus furnishes a method of predicting the value of when 

heat is being transferred by the wall, provided R’ has been previously 
determined from adiabatic tests. 

From Equation (5), it can be seen that physically the recovery factor 
may be thought of as the proportion of the mean stream kinetic energy that 
is recovered at the wall as the air is retarded to zero velocity at the wall 
surface. It is also seen that may be less or greater than Tq depending 

upon whether R’ is greater or less than l.C. 

In general, will not equal since there is a velocity gradient in 

the fluid layers near the walls as the fluid is brought to rest which indicates 
also the presence of a temperature gradient. The temperature gradient in 
the fluid causes heat to flow from the layers near the wall to those further 
removed which are at a lower temperature. The velocity gradient, however, 
causes a transfer of work energy towards the wall in overcoming the viscous 
forces present in the boundary layer. The ratio of heat flow away from the 
wall to the work energy directed toward the wall determines whether T^^ 
will be greater or less than Tg. 

Since the flow of heat away from the wall is controlled by the thermal 
diffusivity = k/y® Cp and the flow of work energy toward the wall is 



controlled by the momentum diffusivity, or kinematic viscosity, 
might be expected that the net effect is controlled by the ratio of ^ to 
which is k//< Cp, the Prandtl number. 

Pohlhausen^"^ in a study of adiabatic wall temperature of an incompress- 
ible fluid flowing over a non-conducting flat plate with a laminar boundary 
layer predicted that the recovery factor was a complex function of only the 
Prandtl number. 

Emmons and Brainerd '' using a differential analyzer to solve the 
differential eqviations of compressible and incompressible flow over a flat 
plate with a laminar boundary layer found that the recovery factor was a 
function of the Prandtl number with some dependence upon the Mach number. 
For Prandtl numbers greater than 0,5 and Mach numbers less than 2.7, 
they predicted that the relationship, R’ = (Pj.)*' would be valid to within 
two percent. 

Eckert and Weise^^ and Ackerman^^ predict that the recovery factor 
is of the order of (Pr) for a turbulent boundary layer. 

Recovery factors for subsonic flow of air in brass tubes were measured 
by Nicolai^^ who reported values of 0,80 to 0.905 for Mach numbers ranging 
from 0.1 to 1,0. Czapek and Marcuse in similar measurements reported 
values of 0.8 5 to 0.96, 

Eber^^ measured surface temperatures on various types of probes 
placed in a supersonic stream and found that for Mach numbers between 
1.5 and 3.0, the recovery factors varied from 0.85 to 0.96 with increasing 



Mach numbers. 



Ketchum, ‘ ^ who compiled the recovery factor data obtained by Jung e 
and Margolskee'^ and Connors and Helfrich^ in previous theses on this 
Project using apparatus similar to that used in this thesis, found that the 
recovery factor for air flowing in a lucite tube at Mach numbers from 
1,0 to 2.6 is a function of the Mach number only, provided the flow is 
turbulent and indicated that the flow became turbulent at inlet diameter 
Reynolds numbers of l.l x 10^. Ketchum’s values of recovery factor, 

0,88 at M = 2.53, 0.84 at M = 1.78 and 0.80 at M = 1.41 were about 3% 
lower than the recovery factors obtained by Eber. At inlet diameter 
Reynolds numbers less than 1.1 x 10^. Ketchum found recovery factors 
that were somewhat erratic but which gave a mean line that for the 
higher Mach numbers agreed with the recovery factors obtained at the 
higher Reynolds numbers, but which at the intermediate Mach numbers 
gave recovery factors two to four percent higher than those obtained at 
the higher Reynolds numbers. At the low Mach numbers the two curves 
remained parallel. 

Friction Factors 

The determination of friction factors in this thesis was undertaken 
to furnish data for the design of any future recovery factor or heat trans- 
fer apparatus and to supplement existing data on friction factors for air 
at supersonic velocities in a tube. 

The apparent friction factor or coefficient is an approximation of 



the true friction coefficient. The latter is defined in terms of the shear 
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stress at the duct wall and is extremely difficult to measure. Previous 
studies indicate the desirability of using the apparent friction factor as a 
convenient form of evaluating results. It is defined in terms of the stream 
pressures, stagnation state and flow per unit area and may be evaluated 
from one -dimensional compressible flow functions for adiabatic flow at 
constant area from the Gas Tables. ^ See Appendix B. 

Experimental results and literature on the subject of apparent friction 

23 

factors in tubes is very limited. Frossel’s" data for supersonic flow 
indicated agreement at fully developed flow with the von Karman-Nikuradse 
relation between friction factor and Reynolds number for incompressible 
flow. Keenan and Neumann^^ also obtain results that show that the friction* 
factor approaches that of incompressible flow at distances in excess of 50 
diameters from pipe inlet. In addition, they show that during the development 
of the boundary layer for a fixed value of L/D greater than 30, the mean 
apparent friction factor decreases as the Reynolds number increases. 

Ketchum, in his analysis of the friction factor data obtained by 
Helfrich and Connors, fotand that the apparent friction factors obtained in 
supersonic flow in textolite and lucite tubes were in good agreement with 
the results obtained by Keenan and Neumann. 



CHAPTER IV 



Investigation Procedure 

The investigation of recovery factors for air flowing adiabatically at 
supersonic velocities resolves itself into a determination of the values 
necessary to the calculation of recovery factors by the following formula: 

~ Taw ~ Tm 3 (^aw/^iig) ~ (^rn/^us) 

Tus - I'm > - (Tm/Tus) 

In order to determine these values the work reported in this thesis 
was started under a plan, the basic features of which were: 

1. The use of dry air at various stagnation temperatures and over a 
wide range of stagnation pressures; 

2. The use of a converging-diverging nozzle to obtain supersonic 
velocities with a maximtim Mach number of about 2.6; 

3. The passage of the air at supersonic speeds through a cylindrical 

tube with about 41 and with facilities for measurement of 

wall temperatures and static pressures at regular intervals along 
its length; and 

4. The reduction to a minimum of heat transfer to the supersonic 
stream by providing thermal insulation in the forms of low con- 
ductivity plastic materials and of an evacuated chamber surrounding 
the lucite cylindrical tube. 

In view of the results obtained by previous investigators whose work 
apparently led to the conclusion that the recovery factor is a function of the 
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stagnation temperature, it was considered desirable to conduct experi- 
mental runs at the following temperature levels: 32°, 65°, 110°, and 130°F. 
In addition, in order to cover a reasonably wide range of Reynolds numbers, 
it was decided to vary the upstream stagnation pressure from about 8 psia 
to 80 psia, the upper limit having been imposed by the inability of a down- 
stream ejector to handle greater flow rates. In all runs critical flow 
through the nozzle was maintained by the use of a steam-driven air 
ejector which maintained very low pressures in the section of the apparatus 
downstream of the nozzle. 

In order to reduce the effect of condensation shock on the static 
pressure the specific humidity of the air entering the test apparatus was 
maintained at about 1.5x10 . It has been estimated that the presence 

of this amount of water vapor would cause a condensation shock v/hich 
woxild affect the static pressures in the negligible amount of about 1 part 
in 100,000.^^ 

A pressure tap was located in each of the eleven bosses spaced at 
two-inch intervals along the test section, and since the pressure gradient 
across the boundary layer can be considered negligible, the pressure 
recorded was essentially the mean stream static value. As described in 
Appendix B these pressure readings (P^,) were used to evaluate the Mach 
Number of the flow at each station, from which the values could 



be determined. 
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Of the three thermocouples located in each boss one was used to 
measure the temperature at a point 0,089 inches radially outward from the 
tube wall and the other two were used to measure the temperature differ- 
ential between two points displaced 0.45 inches radially. From the 
differential value the radial temperature gradient (degrees per inch) was 
determined and the wall temperature was computed by extrapolation. In 
the absence of any heat transfer to the test section, the wall temperature 
thus measured would very closely approximate the adiabatic wall tempera- 
ture at that point. 

Since it was suspected that the previously obtained conclusion 
regarding the variation of recovery factor with the upstream stagnation 
temperature was caused by the fact that the tests were not actually 
adiabatic, it was deemed imperative that the heat transfer be rendered 
negligible. As a result the design and recommendations set forth in 
reference 4 were used in an effort to obtain the predicted temperature 
differential across the bosses of 0.001°F. 

Recovery factor runs were made at the following stagnation tempera- 
tures: 32®, 65®, and 110®F; but preliminary calculations indicated that the 
results wovild be similar to those obtained by Helfrich and Connors, namely, 
that the recovery factor varied with stagnation temperature. It was, 
therefore, found necessary to disregard the twenty-one runs made at 32® 
and 65® and to devote further effort to the reduction of heat transfer. The 
fourteen runs made at 110° were retained and their results submitted for 



I 



IV- 4 



the reasons indicated in Chapter VI. 

Subsequent investigation indicated that the heat transfer to the test 
section could be substantially reduced by effecting the modifications to the 
apparatus as described in Appendix K. Preliminary calculations of 
recovery factors obtained with the modified apparatus indicate that the 
values for a stagnation temperature of llO^F will not be appreciably 
affected but that those at other temperatures will be modified in such a 
direction as to show that the recovery factor is independent of the stagna- 
tion temperature. A more complete discussion of these preliminary 
results can be found in Chapter VI; however, the final results will not be 
available until reported by the investigators who are continuing the work 
of this project. 
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CHAPTF.^ V 
Description of Apparatus 

The apparatus used in this thesis is largely that originally designed by 

0^15 y 

Klingensmith'^ and Wyant and assembled by Junge and Margolskee^ except 

'A 

for the recovery factor test section which was designed mainly by Shoulberg 
and assembled by the authors with some changes. The major improvements 
of the present apparatus over that used by Junge and Margolskee and Connors 
and Helfrich^ for determining recovery factors are: 

(1) Use of lucite test section similar to that used by Connors and 
Helfrich but surrounded by a vacuum chamber to thermally in- 
sulate the test section and to reduce the time required to reach 
thermal equilibrium. 

(2) Use of differential thermocouples to measure temperature 
gradients in the station bosses of the test section in order to 
more accurately determine adiabatic wall temperatures, 

(3) improved method of construction and installation of the test section 
thermocouples, 

(4) Use of a shorter nozzle to reduce the boundary layer effect upon 
determination of recovery factors, 

(5) An improved method of joining the nozzle to the test section to 
reduce misalignments to a minimum, 

(6) Alterations to the pressure measuring system to provide a method 
of measuring absolute pressures at wall of test section by means 



of DC7 03 Silicone oil manometers. 
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(7) Installation of a standard orifice and alterations to the downstream 
piping to provide sufficient straight lengths of uninterrupted piping 
for the orifice. 

(8) Installation of a source of dry air from the Low Temperature 
Laboratory. 

(9) Installation of a copper coil with a larger I.D. and a larger coil 
radius in the constant temperature bath to reduce the pressure 
drop in the coil. 

(10) Installation of several pressure gages and valves throughout the 
system to reduce the pressure losses and to aid in the control of 
the air and the testing of the system for leakage. 

(11) Installation of piping to permit calibration of nozzle and orifice 

by means of gasometer flow. (Piping is uncoupled at the apparatus 
when not in use and does not appear in any of the illustrations of 
this chapter.) 

A brief description of the components of the system appears in the 
following paragraphs with the exception of the orifice installation which is 
described in Appendix J. A recent modification to the test section, made 
after the data for this thesis was obtained, is described in Appendix K. 
Drawings and pictures of the apparatus, as it appeared prior to Run F -36 
and prior to the modification just mentioned, follow the description. 



Air Supply: 



Air is available from four sources; 



Approximate 

Source Maximum Pressure Capacity 

(psiy (CFM) 

Ingersoll-Rand, steam driven, 2 stage 180 124.0 

Chicago , ” ”1 •' 100 69.6 

Joy , motor ” 2 ” 200 250.0 

Low Temperature Laboratory 150 

The last source is obtained from the excess dry air of the air liquefy- 
ing process in the Low Temperature Laboratory which uses the Joy air 
compressor for its air supply. Only air from the Ingersoll-Rand compressor 
or dry air from the Low Temperature Laboratory was used in the r\ms of 
this thesis. 

Air dehumidification: 

Refer to Figure I. Air from the Ingersoll-Rand compressor is passed 
through a coil immersed in an ice -water bath and a filter composed of glass 
wool and copper sponges to femove dirt and to condense oil and water vapor. 
The filter is equipped with a drain to vent the condensate to the atmosphere. 
The air then passes through a trap cooled by the ice-water bath and passes 
upward through a three -tube counter flow, regenerative cooler whose 
temperature approaches -60*^F and in which the remaining water vapor is 
condensed out of the air in a solid phase. Leaving the top of the regenerator, 
the air enters a coil immersed in a dry ice-alcohol bath which cools the air 
several degrees, then passes through a filter of glass wool and copper 
sponges where any solid phase is removed and returns downward through 
the regenerator, A thermocouple (T^,) measures the temperature entering 



200 psi - Joy Compr. 

100 psi - Chicago Corapr. 
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the filter; and to prevent moisture condensation shocks, air from the 
dehumidifier is not admitted to the test section v/hen is higher than 
-40°F. During its downward path, the air removes heat from the upward 
flowing air and emerges from the regenerator as warm dry air. The de- 
humidification process is by -passed when air from the Low Temperature 
Laboratory is used since it has already been dehumidified in the air- 
liquefying process. 

Air Flow and Temperature Control: 

Air from the regenerator or from the Low Temperature Lab now 
passes through the control valves which consist of a needle and a stop 
valve in parallel with a larger globe stop valve. Fine control of the air 
flow is obtained by opening the large valve to the approximate upstream 
pressure desired and then using the needle valve to maintain the exact 
pressure, A by-pass valve is also provided just before the control valves 

r 

to vent excess air to the atmosphere through a muffler. A shell -and -tube 
type Ross air heater, through which a mixture of steam and water is 
circulated, is used to heat the air to within approximately 5®F of the 
temperature of the constant temperature bath through which the air pro- 
ceeds next. This is a large water bath in which the air is heated to approxi- 
mately the desired upstream stagnation temperature as it passes through a 
copper coil (40’ \incoiled length) immersed in the water. The water is heated 
by means of steam to the approximate temperature desired and the tempera- 
ture is then regulated by means of an electric heating coil in the bath with 
an automatic temperature regulator thermostatic control which operates to 
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maintain the desired temperature within 0.2°F. Air is bubbled through the 
bath to stir the water and make the temperature uniform. After leaving the 
bath the air enters the upstream stagnation tank (3” diameter) in which the 
air is slowed to a uniform low velocity by means of baffles and wire screen- 
ing. In the tank are nine thermocouples distributed in a plane perpendicular 
to the flow and equipped with radiation shields to prevent radiation to the 
tank walls. The thermocouples are used to measure the stagnation tempera- 
ture (Tyg) and the temperature distribution of the air in the stagnation tank. 

A pressure tap provides a means of measuring P^g. 

Recovery Factor Apparatus; 

The air from the upstream stagnation tank enters the recovery factor 
apparatus through a convergent textolite reducer (Figure 3) which has an 
exit diameter of 1.25”. The reducer, upstream end plate and the nozzle are 
bolted together firmly with Allen head steel bolts. A stainless steel 
cylindrical ring (Figure 6) with a force fit in the upstream end plate connects 
the reducer to the nozzle. The ring is machined so that there is a near 
perfect transition with no discontinuities to interrupt the flow from the 
reducer to the nozzle. 

The convergent -divergent nozzle (Figure 4) is made of stainless steel, 
with an entrance diameter of 1.25”, a throat diameter of 0.2685”, and an 
exit diameter of 0.5015” and was designed for a Mach number of 2.6 and 
uniform parallel flow at the exit. A flange is provided at each end of the 
nozzle for securing it to the upstream end plate and reducer and for 



fastening the nozzle to the test section. The walls of the nozzle are as 



thin as possible but still provide the necessary strength. Two thermo- 
couple holes are provided to measure temperature distribution along the 
wall but the thermocouples were not installed. 

The nozzle, collar and lucite test section are bolted together firmly. 
The collar is provided to insure nearly perfect alignment of the nozzle 
exit and the test section entrance areas. Special machining and hand 
fitting was required in maniifacturing the lucite collar to make the junction 
between the nozzle and test section as smooth as possible. The recovery 
factor apparatus was designed to be disniantled without the necessity of 
unfastening the nozzle and collar from the test section. Means of re- 
aligning the three parts is provided, however, should such an eventuality 
arise. 

The lucite test section (Figure 5) is similar to that used by Connors 
and Helfrich in regard to diameter (0.5015”), length (21.020”), and loca- 
tion of eleven test stations at 2” intervals along the tube. Lucite was 
used since it is a resin with a low thermal conductivity and reduces the 
effects of longitudinal heat transfer in the test section to a minimum. At 
each station, 1 .88 '‘-diameter bosses are provided to afford a means of 
anchoring the pressure tubing compression fittings to the lucite. In between 
the bosses the tube thickness is reduced to l/8” to further restrict 
longitudinal heat transfer along the tube.. In order to minimize the effects 
of turbulence on pressure measurement at further downstream taps 
the pressure taps are arranged around the tube so that only the first and 
eleventh pressure tap holes are in alignment. Three thermocouple holes 
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are provided at each station. One thermocouple is used to determine the 
wall temperature while the other two, located at different radial distances 
from the tube wall, indicate the temperature gradient in the boss. The 
interior of the lucite tube was made as straight and smooth as possible 
and was given a high polish during manufacture. At the downstream end 
the test section was bolted to the downstream end plate which is also made 
of lucite. 

Surrounding the test section and bolted to the upstream and downstream 
end plates is a thick-walled outer tube (Figure 7) made of textolite for 
strength and low thermal conductivity. In the outer tube are two pipe con- 
nections to the vacuum system, used to evacuate the space between the 
outer tube and the test section, and a smaller pipe connection for the McLeod 
gage vacuum -measuring system. 

The upstream end plate made of lucite contains compression fittings for 
securing the pressure tubing and a means of leading thermocouple wires 
out of the apparatus, (See Figure 6). 

The recovery factor apparatus is sealed throughout with rubber gaskets 
and high vacuum grease at every joint to obtain as high a vacuum as possible. 
Glyptal and air-drying varnish are used externally to seal all joints and 
fittings. 

Downstream air path; 

The air leaves the recovery factor apparatus through a 4*’ I.D. sheet 
steel adapter fastened between the apparatus' and the downstream stagnation 
tank. This tank like the upstream tank contains baffles and wire screens to 
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Figure V-2 
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slow the air to a uniform low velocity. Nine thermocouples and a pressure 
tap are used to determine stagnation temperature (T^jg), temperature dis- 
tribution, and P(jg. Leaving the stagnation tank, the air enters a 2 l/2” I.D. 
pipe leading to the orifice and from thence to the steam driven air ejector 
located in the Steam Laboratory. The air ejector reduces the pressure in 
the test section below atmospheric thus insuring supersonic flow through 
the greater part of the test section even at high rates of flow. 

Vacuum system; 

The vacuum system (Figure 1) used for evacuating the recovery factor 
apparatus consists of an all metal diffusion pump, a trap and a Cenco 
Megavac vacuum pump. The diffusion pump is a fractionating, water 
cooled, electrically heated vacuum pump using octoil and is capable of 
maintaining an ultimate vacuum of 10“^ mm Hg. The Megavac pump 
requires a l/2 hp motor and is capable of maintaining a vacuum of 
0.1 X 10 mm Hg. A ti-ap is provided between the two pumps to prevent 
accidental flow of octoil into the Megavac pump. A flexible coupling is 
provided between the diffusion pump and the apparatus, which effectively 
reduces the vibrations transmitted from the Megavac pump to the apparatus. 
All piping joints are soft soldered and sealed with glyptal and air-drying 
varnish with the exception of two unions which are sealed with high vacuum 
grease and glyptal* The best vacuum obtained was 1.81 X 10“-^ mm Hg. 

Pres sure -measu ring system; 

This system provides mercury and silicone oil (Dow-Corning 7 03) 
manometers as shown in Figure 10. Provided the pressure in the absolute 
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manometer manifold is less than C.l mm Hg, absolute pressures may be 
read directly from the silicone manometers which have one leg connected 
to the absolute manometer manifold which is evacuated by means of a 
Cenco Hyvac Pump capable of maintaining 3 x 10“^ mm Hg pressure. All 
silicone manometers are equipped with stop cocks to prevent oil from 
entering the pressure lines or the absolute manometer manifold. One 
manometer containing silicone oil is provided to compare the pressure 
in the absolute manometer manifold with that in the recovery factor 
apparatus. A McLeod gage calibrated by Jung e and Margolskee is used 
to determine the pressure in the recovery factor apparatus. A Bourdon 
pressure gage is provided for measuring P^s than 43 psia. 

Silicone oil is used because of its low vapor pressure, low air solubility, 
and a specific gravity (1.09) near that of water. It also does not dissolve 
the grease used to keep the stop cocks air tight. 

Temperature-measuring system: 

The upstream and downstream stagnation tank thermocouples are 
those installed and calibrated by Jung e and Margolskee. The thermo- 
couples installed in the test section consist of eleven used to measure sta- 
tion wall temperatures and calibrated as in Appendix H and twenty -two 
used to measure the temperature gradient across each station boss. The 
thermocouples were made of #30 constantan and copper wire twisted 
together for a length of l/8” and soft soldered. The junction was care- 
fully placed in the drilled hole of the station boss and the hole filled with 
high vacuum grease using a hypodermic needle. The 44 wires were led 



out of the apparatus through 45° holes in the upstream end plate and thence 
between the end plate and the textolite outer tube. An effective seal was 
obtained by means of high vacuum grease, rubber gaskets, and externally 
applied glyptal and air -drying varnish. The cold junctions of the station 
wall thermocouples were placed in an ice water mixture contained in 
thermos bottles along with the cold junctions of the stagnation tank thermo- 
couples. The temperature reading was obtained by means of a sensitive 
galvanometer and potentiometer in conjunction with a calibrated standard 
cell and a 6 volt storage battery. A switchboard was provided for selection 
of the desired thermocouple reading. 
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CHAPTER VI 
Discussion of Results 

General 

The results presented here were obtained from fourteen experimental 
runs, F-1 to F -9 and F-17 to F-21, using the apparatus described in 
Chapter V. Each run was made with the upstream stagnation temperature 
and pressure held constant while measurements were made at eleven sta- 
tions along the test section to obtain data necessary to the determination of 
recovery factors, friction factors, Mach numbers, and Reynolds numbers. 
All r\ins were conducted at an upstream stagnation temperature of about 
110°F and at upstream pressures varying from 8 to 80 psia, the range of 
inlet diameter Reynolds numbers being 0.47 x 10^ to 4.8 x 10^. 

Rtins F-17 to F-21 were carried out as duplicating runs to insure 
that the results obtained in Runs F-1 to F -9 could be repeated. In general, 
there was excellent agreement between the data obtained from runs which 
duplicated each other. 

The results of the experimental work are presented in tabiilar form 
at the end of this chapter and are shown graphically in Figures VI -1 to 
VI-10. Except for Figure VI-2 only data for points upstream of a 
suspected shock were plotted. 

Figures VI-! and VI-2 show the distribution of wall pressures and 
calculated adiabatic wall temperatures along the test section for the 
first nine runs. From these plots can be seen the influence of shocks 
near the downstream end of the cylindrical tube as well as the effects of 
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heat transfer to the test section at its inlet and outlet ends. 

In the pressure-ratio plot, Run F-9, which was made with the lowest 
inlet diameter Reynolds number, shows a pressure ratio which increases 
along the length of the tube, a situation which is impossible in the one- 
dimensional theory used as a basis for the work of this thesis. Since 
similar results were reported by Kaye, Keenan, and McAdams and 
Ketchum,^^ it may be concluded that the discrepancy was caused by the 
limitations of the theory rather than by errors in measurement. 

In the temperature -ratio plot Run F-5 shows a distribution which is 
inconsistent with that of other similar runs. In checking through the 
original data used to calculate the results of this run, it was found that 
this was the first of five runs carried out in one day. As a result it is 
believed that temperature equilibriiim had not actually been attained 
when readings were commenced and that, consequently, all results of 
Run F-5 should be disregarded. 

Recovery Factors 

The results which are of primary importance in the determination 
of recovery factors are those relating the recovery factor to the Mach 
number; they are shown in Figure VI-3, which includes the results of 
runs made over the entire range of inlet diameter Reynolds numbers from 
0.47 X 10^ to 4.8 X 10^. Discounting the values for Run F-5, the deviation 
of any of the points from a mean line passed through them is less than 
one percent of the recovery factor. 
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In attempting to correlate the results presented here with those of 
previous investigators, a comparison with the data obtained by Ketchum 
was made. His data led to the apparent conclusion that the recovery 
factor depended on the inlet diameter Reynolds number as well as the 
Mach number; as a result he divided his results into two regimes, the 
line of demarcation being an inlet diameter Reynolds number of 
1.05 X 10^. Ketchum concluded that only in the region of higher inlet 
Reynolds could the recovery factor be regarded as a function of Mach 
number only. The results presented here do not indicate the necessity 
for such a division since the recovery factors determined seem to be 
independent of the Reynolds number. In addition, there is considerably 
less scatter in the region of low flow rates. 

At low values of Mach number the results indicated in Figure VI-3 
are in excellent agreement with those obtained by Ketchum in the regime 
of high rates of flow, but at the higher Mach numbers there is a signifi- 
cant deviation with Ketchum’s recovery factors being somewhat lower. 
The net result of this difference is that the mean of the values presented 
here is more nearly linear. 

In Figure VI-4 recovery factors are plotted against diameter 
Reynolds number for each run. The plot is interesting in that the lines 
connecting successive stations for each run are roughly parallel at a 
slope of 0,47. This suggests a strong dependence of recovery factor 
on viscosity and, consequently, on the Prandtl number; however, it has 
been impossible to correlate the results indicated in Figure VI-3 with 



simple formiilas involving fractional exponents of the Prandtl number. 

The plot of recovery factor vs. length Reynolds number (Figure VI-5) 
shows two flow regimes which can best be characterized by the inlet 
diameter Reynolds number. If the results of Run F -9 are disregarded on 
the basis that the one -dimensional theory is inadequate for the conditions 
of the run, then one regime could be identified with flow for which the 
inlet Reynolds number is between 0.83 x 10^ and 1.8 x 10^ and for which 
the recovery factor could be expressed as a function of length Reynolds 
number. In the other regime for inlet Reynolds number greater than 
1.8 X 10“ the recovery factor seems to be independent of the length 
Reynolds number and to be more dependent on the L/D ratio. 

Friction Factors 

Although local apparent friction factors were calculated along with 
the mean apparent values, the preliminary plots of the former showed 
no consistency and as a result, only the latter values have been plotted 
and will be considered in this discussion. The local apparent values, 
however, are included in the tables of results. 

The most significant results obtained are those indicated in 
Figure VI -6, which is a plot of friction factor as a function of length 
Reynolds number. The plot is quite analogous to the corresponding 
relationship between friction factor and diameter Reynolds number for 
incompressible flow. It clearly defines three flow regimes which appear 
to be essentially laminar, transitional, and turbulent. In this figure the 
region of transition lies roughly in the range of Reynolds numbers from 
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15 X 10^ to 3u X 10^; this compares reasonably well with Ketchum’s data 
which indicated that the limits of this regime were 7 x 10^ to 25 x 10^. 

The friction factor reaches a minimum of about 0.0013 as it reaches 
the transition region and after transition decreases steadily with increasing 
Reynolds number from a maximum value of about 0.0035, 

In the plot of friction factor vs. L/D ratio (Figure VI-7) the results 

1 7 

shown are in good agreement with the results of Keenan and Neumann 
and Ketchum except for the values at the second station (l/D of about 5). 

In this report the frction factors for runs conducted at the higher inlet 
Reynolds numbers increase in going from station 2 to station 3 whereas 
in the other reports the friction factor decreases. Since the friction 
factor is dependent only on the Mach number, or the pressure ratio, it is 
not immediately evident why the difference exists. 

Figure VI -8, which is plotted to a linear scale, again shows the 
presence of three flow regimes with a transition occurring at a diameter 
Reynolds number of about 10^. After transition the friction factor appears 
to level off at a value of about 0.003 or decrease slightly from that value 
as the diameter Reynolds number increases. In addition, in the region 
beyond transition, the friction factor seems to be reasonably independent 
of the L/D ratio. 

Effect of Va cuum System 

Although no data was observed for the exclusive purpose of deter- 
mining the effect of the evacuated chamber on heat transfer, it is 
interesting to consider the apparent results obtained. Runs F-2 and F-18 
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were conducted at upstream stagnation pressures which differed by less 
than two percent and at upstream stagnation temperatures which differed 
by less than 0.1 percent, but F -2 had a of 0.0035 mm Kg and F-18 of 

0.0091 mm Hg. In this case the recovery factors of F-18 v/ere roughly 
one-half percent higher than those of F-2. Unfortunately no other set of 
duplicating runs was made under such conditions as to afford a comparison 
as good as the above since in most cases the temperatures and pressures 
differed by greater amounts while the values of Pvac showed a smaller 
variation. 

Preliminary Results with Modified Apparatus 

After the recovery factor apparatus was modified to reduce heat 
transfer to the test section (See Appendix K), runs were continued in an 
effort to investigate the effect of stagnation temperature on recovery 
factor. The initial runs were made at stagnation temperatures of about 
110°F to check the validity of the results reported here. In Figure VI-10 
the results of two of these runs are compared with data obtained from 
the original apparatus for approximately the same diameter Reynolds 
nximber at inlet. It can be seen that the recovery factor for the more 
recent runs is lower by about one percent at the higher Mach numbers 
and that at lower Mach numbers the difference is roughly one -half percent. 

The data obtained with the modified apparatus, consequently, more 
closely approximates that obtained by Ketchum despite the fact that his 
apparatus was more susceptible to heat transfer than that used for this 



report. However, since the differences are of such a low order and since 
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the adiabatic v/all temperature for stagnation temperatures of 110°F should 
be sufficiently close to room temperature to render heat transfer effect 
negligible, it is felt that the results presented here are sufficiently valid 
for the purposes for which they will be used. 
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Rur F-36 



Pus “ 24,91 psia 
Tus - 569.57 ^Fabs 
Pvac ~ 0.011 mm Hg 
Rcd at inlet - 1.508 x 10^ 



STN Pw/Pus ^ Taw/ Tus 



1 


.03970 


2.640 


.9342 


.8870 


2 


.04360 


2.488 


.9327 


.8782 


3 


.04589 


2.406 


.9353 


.8796 


4 


.05287 


2.164 


.9350 


.8655 


5 


.06291 


1.945 


.9370 


.8537 


6 


.07138 


1.778 


.9391 


.8428 


7 


.08146 


1.614 


.9424 


.6318 


8 


.12935 


1.121 


.9501 


.7516 



Run F -43 

Pus " 24.28 psia 
Tus ” 565.02 °Fabs 
Pyaj. - 0.0046 mm Hg 
Rej^ at inlet - 1.457 x 10^ 

STN Pw/Pus 


M 


Taw/ Tus 


R’ 


1 


.04088 


2.591 


.9349 


.6864 


2 


.04490 


2.440 


.9314 


.8738 


3 


.04745 


2.353 


.9349 


.8761 


4 


.05383 


2.164 


.9346 


.8652 


5 


.06314 


1.940 


.9368 


.8528 


6 


.07 084 


1.788 


.9392 


.8442 


7 


.08259 


1.597 


.9427 


.8304 


8 


.13151 


1.105 


.9511 


.7 507 



CI-iAFTEP VII 



Conclusions 

For air flowing in a tube at supersonic velocities the recovery factor 
is very nearly a linear function only of the Mach nunober in the range of 
Mach numbers from 1.3 to 2.6 and inlet diameter Reynolds numbers of 
0.47 to 4.8 X 10^. The recovery factor varies from 0.81 to 0.905 at Mach 
numbers of 1.3 and 2.6 respectively. 

Preliminary results with the apparatus modified to reduce heat 
transfer effects to a minimum indicate that the recovery factor is inde- 
pendent of stagnation temperature. In addition the results indicate that 
the recovery factors measured with the modified apparatus are about 
0.5% lower than those reported in this thesis. 

Results obtained in this thesis show that it is possible to define 
three flow regimes in terms of length Reynolds number with a transition 
regime occurring between the Reynolds numbers of 15 x 10^ and 30 x 10^. 

In supersonic flow in a tube, the results show that variation of the 
mean apparent friction factors with length Reynolds number is analogous 
to the relation between friction factor and diameter Reynolds number in 
incompressible flow. The mean apparent friction factor reaches a 
minimum value of 0.0013 as it reaches the transition region and after 
the transition region it decreases steadily with increasing Reynolds 



numbers from a maximum value of about 0.0035. 



CHAPTER VIII 



Recommendations 

Based on the results and experience obtained in preparing this thesis, 
the follov/ing recommendations are made: 

1. Further runs at stagnation temperatures other than 110°F should 
be made to determine whether the recovery factor is a fvinction of the 
stagnation temperature, 

2. The range of Mach numbers covered should be extended by using 
a supersonic nozzle designed for a Mach number of about 4.0. 

3. In order to increase the range of flow rates which can be 
accommodated by the apparatus, the capacities of the ejector and re- 
generator should be increased. 

4. The actual effect on heat transfer of the evacuated chamber 
surrounding the test section should be investigated to determine whether 
vacuums of the order obtained in this thesis are sufficient to render heat 
transfer negligible. 

5. In order to eliminate the effect of heat transfer to the air flowing 
between the stagnation tanks, it is suggested that a study be made of the 
possibilities of enclosing in an evacuated chamber the apparatus from the 
inlet of the upstream stagnation tank to the outlet of the downstream 
stagnation tank. It is further suggested that a means be provided for 
temperature control of the external surface of the test section, 

6. A.n attempt should be made to correlate the experimental data 
by using Reynolds numbers based on a viscosity which is evaluated at 








m 









some temperature between the mean stream and the wall temperatures. 

7. For measuring upstream stagnation pressures greater than 
43 psig, a gage which is more accurate than the Bourdon gage now in 
use should be installed. 

8. The supersonic nozzle used in this investigation should be re- 
calibrated while it is not assembled in the recovery factor apparatus in 
order to better cover the range of expected nozzle throat Reynolds numbers. 

9. A filter should be installed in the pipe line through which passes 
the dry air supplied by the Low Temperature Laboratory. 

10. Safety traps should be installed in the low pressure legs of the 
absolute manometers in order to prevent the loss of silicone oil through 
flow into the vacuxxm system manifold. 

11. The constant -temperature bath should be provided with a heater 
of greater power output since the one currently in use is incapable of 
maintaining a constant temperature of 110°F. 

12. A stop valve should be installed in the section of piping 
immediately upstream of the Cenco-Megavac pump to facilitate the 
detection of vacuum system leaks through the pump. 

13. The use of solid nylon as a material for future test sections 
should be investigated. Mr. Chase, of the Cummings Machine Works, 
Boston, has shown that nylon in this form has a greater machinability 



than lucite. 



APPENDIX A 



Brief History of Project DIC 6418 



The work reported in this thesis was carried out as a part of 
Project DIC 6418, which is sponsored by the United States Navy through 
the Office of Naval Research, Since the Project has been in existence 
over a period of several years and since the results of previous investi- 
gators are directly related to those now being reported, a brief history 
is presented here. 

The following is quoted from reference 4, the publication of which 
immediately preceded this phase of the Project: 



“The basic apparatus used to measure both recovery factors and 
heat transfer coefficients was designed in 1946 by Klingensmith (5), His 
work included the design of a supersonic (M = 2,5) nozzle, the design of 
a recovery factor test tube made of textolite to which the nozzle would be 
directly attached, and some calculations for the design of a heat transfer 
test section. The heat transfer test section was completely designed by 
Wyant (3) in 1947. 

“Junge and Margolskee (6) began construction of the Klingensmith 
apparatus in February, 1947. The supersonic nozzle and textolite tube 
were machined at the Boston Navy Yard, The nozzle was calibrated by 
Ketchum (7) (who was in charge of experimental work between February 
1946 and November 1948) in early March, 1947. In early July, 1947, 
tests were made on the recovery factor apparatus. Prior to the assembly 
of the apparatus small burrs were visible in the test section at stations 
3 and 7, When test runs were made, the pressure distribution was not 
smooth, so the inside of the test tube was burnished with a glass rod to 
smooth the edges of the pressure taps. At the same time the test section 
was shortened by machining off two inches of the tube in order to give a 
greater length of supersonic flow. Further tests showed that station 5 
still gave an vinsatisfactory pressure reading, so the hole was sealed with 
glyptal and a new pressure tap was drilled l/8 in, further downstream. 

At the same time a brass collar was machined to fit between the nozzle 
and the test section. Additional trial runs were made during which it 
was observed that the pressure distribution was satisfactory near 
station 5 when Pus was 21 psia, but unsatisfactory when Pus was 96 psia. 



The nozzle was then rotated 90*^ relative to the tube, but no appreciable 
change in pressure distribution was observed. Recovery factor runs 
JMl-9 were then made (6). 

*Tn October, 1947, Connors and Kelfrich (4) took over operation of 
the apparatus and made runs T 1 to T7 to verify the results of Jung e and 
Margolskee (6). A lucite tube was machined and then installed. This tube 
had eleven stations and an extension v/hich could be installed if desired. 

A diffuser v/as also provided, since it was thought that the diffuser v/ould 
increase the length of supersonic flov/ in the tube. The diffuser did not 
perform satisfactorily except at very low upstream stagnation pressures 
(belov,^ 100 cm Hg absolute). At about the same time the manometer 
system was revised. 

“Runs LI to LZO were then made (4),,v/ith a range of stagnation 
temperatures between 33° F and 140° F. It appeared that the recovery 
factor might be a function of the upstream stagnation temperature, but 
careful investigation showed that the apparatus used by Connors and 
Helfrich was reliable only when t^s was about 105° F, so no conclusions 
could be drawn about a Tus - R’ relationship, 

“In June, 1948, the heat transfer section was installed in place of 
the lucite recovery factor section by Larkin, Shoulberg, and Ketchum, 

At the same time the design of an improved recovery factor apparatus 
was begun. ’’ 

During the period from June 1948 to May 1949 Shoulberg, Larkin, 
and England conducted an investigation of heat transfer coefficients for 
air flowing in a tube at supersonic velocities. Runs K-1 to K-28 were 
made with the apparatus designed by Wyant. Runs S-IA to S-3B v/ere 
made after the apparatus had been modified to permit variation of the 
wall temperature without changing the air flow conditions. The super- 
sonic region investigated lay in the range of Mach numbers from 
1.0 to 2.5 and diameter Reynolds numbers from 30,000 to 400,000. 

In addition, Shoulberg, Larkin and England completed the design 
of an improved recovery factor apparatus which was subsequently 
modified and used for the experiments of this thesis. 
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Working at the same time as the above, Ketchum, in a thesis for the 
doctorate degree, compiled the recovery factor and friction factor data 
obtained by Junge and Margolskee and Connors and Helfrich with the heat 
transfer coefficients obtained by Shoulberg and made a more thorough 
analysis of the data obtained, Ketchum also made a very complete review 
of the literature pertaining to this field. 

The work accomplished on the project from the period from June 1949 
to May 1950 by V/elch, Volonte, Toong, and Young may be summarized 
chronologically as follows: 

1. Design of Foelsch and Sauer supersonic nozzles. 

2. Revision of drawings and awarding of contract for the manufacture 
of the recovery factor apparatus. 

3. Making of subsonic heat transfer runs followed by dismantling of 
heat transfer apparatus. 

4. Assembly of new recovery factor apparatus without the test section 
thermocouples. 

5. Installation of pressure gages and valves in the air piping. 

6. Overhauling of regenerator and renewal of all filtering material, 

7. Design and installation of an orifice plate. 

8. Calibration of nozzle and orifice plate with the gasometer. 

9. Manufacture, calibration, and installation of test section 
thermocouples and wiring. 

10. Installation of larger copper coil in constant temperature bath. 



11, Design and assembly of vacuum system for the recovery factor 
apparatus and modification of the pressure measuring system. 
IZ. Making of recovery and friction factor runs F-1 to F-35. 

13, Modification of the recovery factor apparatus to reduce amoxznt 
of heat transfer between room and test section. 

14. Making of experimental runs F-36 to F-45 to determine the 



effects of the modification. 



APPENDIX B 



Calculation of Recovery Factors and Friction Factors 
Recovery Factors 

In general the method of calculation of recovery factors follows that of 
Helfrich and Connors.^ By definition the recovery factor for a perfect gas 
is given by: 

j^* _ ^aw " r 1 - ~ Taw 

^us “ ^us " ^m 

Since T^g is directly measurable and can be closely estimated, it 

remains only to determine in order to compute the recovery factor. 

For a gas flowing isentropically to a supersonic nozzle and adiabatically 
through the nozzle and duct of constant area, it can be shown that the product 
(Pw/^us) (A/A*) is the same as for the case of isentropic expansion. The 
pressure ratio, Mach number, and flow coefficient are related by the 
expression^ 




Since (A/A*) is fixed in this case, the Mach number can be determined when 
the flow coefficient and pressure ratio are known. Using tabulated values 
of (Pw/Pus) (a/a*) as a function of Mach number from the Gas Tables, ^ 
a plot of (l/Cw) (Pw/Pus) vs, M was constructed and is shown in Figure C-1, 
The pressure ratio at a given station can be evaluated from the observed 
data, and the flow coefficient can be determined from the nozzle calibration 





«i 



t 



B-2 



plot of vs Rej 3 * (Figure C-2) after Rej^* has been connputed from the 
following: 



Rej)* = 



(w/A'!=) D’f' 



where D=<‘= is a constant. 



From the isentropic tables for k = 1,4 

tVTus = 0.83333. 

Having computed T* the viscosity at the nozzle throat can be determined from 
the temperature-viscosity plot (Figure C-3). From the equation for maxi- 
mum isentropic mass rate of flow for a perfect gas through a nozzle: 

(ws/A*) = 0.03012 

sj 

the isentropic flow rate is readily determined. In calculating the throat 
Reynolds number for the purpose of determining the flow coefficient, it is 
sufficient to use the isentropic rather than the actual value of flow rate. In 
an extreme case where a known flow coefficient of 0.9600 is neglected, the 
error in Rej^* is 5.8% with a resulting error of 0.09% in estimated flow co- 
efficient, which error is well within the limits of precision with which the 
nozzle was calibrated. 

Once the Mach number has been determined from Figure C-1, the mean 
stream temperature can be evaluated from the following expression obtained 
from a one-dimensional analysis of adiabatic flow® 

^ _J 

The calculations are simplified by using the tabulation of the temperature 
ratio as a function of Mach number in the one -dimensional isentropic tables 



of the Gas Tables. 



This method of recovery factor calculation differs from that of Helfrich 



and Connors in that: 

1. The flow coefficient has been included in the expression relating 
pressure ratio, area ratio, and Mach number; and 

2. In accordance with the assvimption of a perfect gas the Joule-Thomson 
effect has been neglected in computing the stagnation temperature at each 
station. 

The recovery factors for Run F -6 have been calculated by the methods of 
Helfrich and Connors and of Junge and Margolskee;^ the results are presented 
for comparison in Figure B-1. 

Friction Factors 

From a one -dimensional analysis of compressible adiabatic flow in a 

Q 

tube of constant diameter it can be shown that 



Since the 0 ftinction is tabulated as a function of Mach number for such flow 
in the Gas Tables, the calculation of friction factors is reduced to a deter- 
mination of the Mach number, which is evaluated incidental to the calctilation 
of recovery factors. 

To compute the local apparent friction factor, i.e., the value which is a 
mean for two adjacent stations of the test section, it is necessary only to 
know the Mach numbers at the two stations. Then 



4 "f 



D 



I - ^ K I 

K 





0 



D 



Letting f be a value which is the mean of £j and 



o 

But ^ max2 distance between any two adjacent stations and 



is constant; therefore 

D 


aL = <t>. - 




^ ( (t). - 

4 



To compute the mean apparent friction factor, i.e., the value which is a 
mean between the first station and any other station, the above expression 
can be transformed to 

f { <t>. - 0.) 

4 

where a. is the distance between the stations involved. 
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APPENDIX C 

Calculation Data, Constants, and Coefficients 
Calculation Data 

Although the data used in the calculation of the recovery and friction 
factors was compiled in tabular form in order to insure consistency among 
the several persons making the calculations, they are appended here in the 
form of graphs. The necessary data include the following: 

1. M vs (1/C^) (P^/P^s). the source of which is described in 
Appendix B and which is shown in Figure C-1. 

2. Cw vs ReD** the source of which is described in Appendix G and 
which is shown in Figure C-2. 

t 

3. yCA. VS T, which data is the same as that used in references 1, 2, 4, 
and 10 and which is shown in Figure C-3. 

4. Pressure conversion data plotted in Figure C-4. 

5. Thermocouple calibration data, the basis for which is described in 
Appendix H and which is tabulated in Table C-1. 

6. '^rr/'^us 'vhich was obtained from the tabulation of T/Tq in 

Table 30 of the Gas Tables (One-Dimensional Isentropic Compressible-Flow 
Functions for k = 1.4). 

vs M, which was obtained from Table 42 of the Gas Tables 
(Fanno Line - One-Dimensional Compressible-Flow Functions for Adiabatic 
Flow at Constant Area with Friction for k = 1.4). 
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Calculation Constants 

The constants used in the calculation of recovery and friction factors 

are; 



A 


- 0.1978 in^ 


R 


- 53.345 ft-lb/lb-°F 


A* 


- 0.05662 in^ 


^s 


- 0.340 in. 


D 


- 0.5018 in. 


^^s 


- 0.450 in. 




- 0.2685 in. 




- 0.251 in. 


g 


- 32.174 ft/sec^ 


^1 


- 2.07 6 in. 


k 


- 1.400 


^2 


- 4.076 in. 


Ll 


- 0.810 in. 


X 3 


- 6.076 in. 




- 2.810 in. 


X 4 


- 8.076 in. 




- 4.81 0 in. 


^5 


- 1 0.076 in. 


1-4 


- 6.8 1 0 in. 


^6 


- 12.076 in. 


L 5 


- 8.8 1 0 in. 


^7 


- 14.076 in. 


1^6 


- 1 0.81 0 in. 


^8 


- 16.078 in. 




- 12.810 in. 


X 9 


- 18.076 in. 


1-8 


- 14.810 in. 


^10 


- 20.076 in. 


^9 


- 16.81 0 in. 


xii 


- 22.076 in. 


^10 


- 13.810 in. 


°Fabs = + 459.69° 


Lfi 


- 20.81 0 in. 







Calculation Coefficients 



Flow Rate 



The maximum isentropic mass rate of flow of a perfect gas is given by 

A* 












R 









Then 



and 



w 



A 



^ = 0 5^2. 



Ws = OOlbOXZ. 



yr 






where Wg is in lb/ sec, p^g psia, and T^g in ®Fabs. 

At the nozzle throat the isentropic flow rate per vinit area is 



w. 



= -A* = 2543.4 Wg. 

where Gg* has the units of lb/sec-ft2. At any station along the test section 

^w 

G = — 728.13 C\y Wg 



Diameter Reynolds Number 

The Reynolds nunriber based on the diameter is given by 



At the nozzle throat 



Rbd = 



GD 



Rej5 



* = 



0.02337 G* 



and at any station along the test section 



Rcd 



0.041 82 G 



C-4 



where the G’s have the units of lb/ sec-ft^ and the viscosities which are 
deternnined at T* and T^^j, respectively, are in Ib/sec-ft. 

Adiabatic Wall Temperature 

The coefficient used in calculating the extrapolated adiabatic wall 
temperature is derived from the following: 

*s “ *aw ^s ~ ^w 

Cits A rs * 

For the constants of this report 



ts - taw 



0.1978 Ats . 




t 



Table C-I 

Thermocouple Calibration Data* * 



mv 






-Lo 




-16.415^ 


-0.9 




-11.415 


-0.8 




- 6.451 


-0.7 




- 1.523 


-0.6 


3.58°F 


3.369 


-0.5 


8.40 


8.225 


-0.4 


13.18 


13.045 


-0.3 


17.94 


17.828 


-0.2 


22.65 


22.576 


-0.1 


27.34 


27.306 


0 


32.00 


32.000 


0.1 


36.64 


36.648 


0,2 


41.25 


41.260 


0.3 


45.83 


45.854 


0.4 


50.40 


50.431 


0,5 


54.94 


54.97 1 


0.6 


59.45 


59.493 


0.7 


63.94 


63.997 


0.8 


68.41 


68.465 


0.9 


72.86 


72.915 


1.0 


77.28 


77.347 


1.1 


81.67 


81.744 


1.2 


86.03 


86.104 


1.3 


90.37 


90.446 


1.4 


94.68 


94.770 


1.5 


98.99 


99.076 


1.6 


103.28 


103.364 


1.7 


107.54 


107.635 


1.8 


111.77 


111.869 


1.9 


115.98 


116.085 


2.0 


120.17 


120.283 


2.1 


124.35 


124.463 


2.2 


128.51 


128.625 


2.3 


132.65 


132.769 


2.4 


136.77 


136.895 


2.5 


140.90 


141.003 


2.6 


144.99 


145.112 


2.7 


149.06 


149.184 


2.8 


153.12 


153.238 


2.9 


157.16 


157.274 


3.0 


161.19 


161.292 



* *See Appendix H, 
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Figure C-3 
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Figure C-4 



APPENDIX D 



Sample Calculation 

Following is a tabulation of the operations involved in the calculation 
of recovery and friction factors. A sample data and calculation sheet for 
Run F -6 is included at the end of this appendix. The constants and coef- 
ficients used in the calculations are evaluated in Appendix C. 



Column Item Dimensions 



(l) L/D none 



(2) x/D none 



(3) 


T room 


°F 


(4) 


^bar 


in. Hg 


(5) 


^bar 


psi 


(6) 


Pus 


cm. Hg 



Operation 

L/D values are constant; the 
values for stations 1-11 are, 
respectively, 1.614, 5.600, 9.586, 
13.571, 17.557, 21.542, 25.528, 
29.514, 33.499, 37.485, 41.471. 
x/D values are constant; the 
values for stations 1-11 are, 
respectively, 4.137, 8.123, 

12.108, 16.094, 20.080, 24.065, 
28.051, 32.037, 36.022, 40.008, 
43.994. 

Observed data. 

Observed data (corrected to 32°F). 
Multiply (4) by 0.49119. 

Observed data (manometer left 
column). 



— o 



Column 


Item 


Dimensions 


Operation 


(7) 


Pus 


cm. Hg 


Observed data (manometer right 
column). 


(S) 


» » 


9 9 


(7) - (6), 


(9) 


1 9 


psig 


Multiply (8) by conversion factor 
for mercury evaluated at 
temperature (3); Figure C-4. 


(10) 


9 t 


psia 


(9) + (5). 


(H) 


Pw 


cm. Hg 


Observed data (manometer left 
column). 


(12) 


1 9 


9 9 


Observed data (manometer right 
column). 


(13) 


9 9 


9 9 


(12) - (11), 


(14) 


9 9 


psig 


Multiply (13) by conversion 
factor for mercury evaluated at 
temperature (3); Figure C-4, 


(15) 


9 9 


psia 


(5) -(14). 


(16) 


9 9 


cm. Silicone 


Observed data (manometer left 
column). 


(17) 


9 9 


9 9 


Observed data (manometer right 
column). 


(18) 


9 9 


9 9 


(17) -(16). 
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Column 



Item 



Dimensions 



Operation 



(19) 


Pw 


psia 


Multiply (18) by conversion 
factor for silicone oil 
evaluated at temperature (3); 
Figure C-4. 


(20) 


^US 


mv 


Observed data; three values 
measured at different times. 


(21) 


t « 


Op 


Convert average of three values 
in (20) using thermocouple 
calibration data. 


(22) 


T 

us 


°Fabs 


(21) + 459.69. 


(23) 


T=«' 


°Fabs 


0.83333 X (22). 


(24) 




lb/ sec-ftx 1 0^ 


Find value of viscosity at 

f 

temperature (23); Figure C-3. 


(25) 


\/'^US 


(OFabs)"^ 


v/(22) 


(26) 




lb/ sec 


0.03012 X (10)/(25). 


(27) 




Ib/sec-ft^ 


2543.4 X (26). 


(28) 


Rejj- 


none 


0.02337 X (27)/(24). 


(29) 




none 


Flow coefficient evaluated at 
(28) from Figure C-2. 


(30) 


^w/^us 


none 


(15)/(10) for mercury manometers 



(19)/(10) for silicone oil 
manometers. 

(31) (l/C^)(Pw/Pus) none (30)/(29). 



JC 




Column 



Item 



Dimensions 



Operation 



(32) 


M 


none 


Evaluated at (31) from Figure C- 


(33) 


Tm/Tus 


none 


Evaluated at (32) from Table 30 
of Gas Tables. 


(34) 


Tm 


OFabs 


(33) X (22). 


(35) 


ts 


mv 


Observed data. 


(36) 




Op 


Convert (35) using thermocouple 
calibration data. 


(37) 


^ tg 


mv 


Observed data. 


(38) 


(t + At)g 


mv 


(35) + (37). 


(39) 


(t + A t)s 


Op 


Convert (38) using thermocouple 
calibration data. 


(40) 


^ tg 


op 


(39) -(36). 


(41) 


” ^aw 


Op 


0.1978 X (40). 


(42) 


*aw 


Op 


(36) -(41). 


(43) 


^us ~ ^aw 


op 


(21) -(42). 


(44) 


T - T 


oPabs 


(22) -(34). 


(45) 


1 - R* 


none 


(43)/(44). 


(46) 


R’ 


none 


1 - (45). 


(47) 


aw/ ^ us 


none 


(42) + 459.69 

JJzf 


(48) 


G 


lb/ sec-ft^ 


728.13 X (26) X (29). 


(49) 




lb/ sec -ft X 1 07 


Find value of viscosity at 
temperature (34); Figure C-3. 


(50) 


Rei5 


o 

1 


0.04182 X (48)/(49). 



Colvimn 



Item 



Dimensions 



Operation 



(51) 


ReD 


X 10-5 


(52) 




X 10-5 


(53) 




X 10”^ 


(54) 


^«x 


X 10"^ 


(55) 




none 


(56) 


is ^ 


none 



(57) f X 10^ 

(58) none 

(59) D/4 X 10^ 



(60) 


f 


X 10 


(61) 


hi 


cm 


(62) 


^2 


cm 



Average of values for adjacent 
stations in (50), 

(1) X (50). 

Average of values for adjacent 
stations in (52), 

(2) X (50). 

Evaluated at (32) from Table 42 
of Gas Tables. 

In (55) subtract value for 
station 2 from value for 
station 1, station 3 from 
station 2, etc. 

0.062725 X (56). 

In (55) subtract value for each 
station from value for station 1 . 
Values are constant; for stations 
1-11 they are, respectively, 
infinity, 6.2725, 3.1363, 2.0908, 
1.5681, 1.2545, 1.0454, 0.8961, 
0.7841, 0.6969, 0.6273. 

(58) X (59). 

Observed data. 



Observed data. 




t 



Column 



Item 



Dimensions 



Operation 



(63) 


hi -h2 - 1 


cm 


(61) - 


(62) - 1. 


(64) 


80.25 - h2 


cm 


80,25 


- ^2* 


(65) 


p 

vac 


mm Hg, 


(63) X 


(64)/975. 
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APPENDIX E 



Design of Nozzle 

In the design of the improved recovery factor apparatus by Shoulberg, 
the nozzle contour designed by Klingensmith and used in the original 
recovery factor apparatus was retained with changes only in the dimension 
reference plane and in the mounting arrangements. 

In an effort to obtain a nozzle of shorter length than the Klingensmith 
nozzle and thus reduce the boundary layer effect upon flow through the noz- 
zle, two nozzles were designed in June- July 1949 by G. J. Van Wylen, 
research assistant to the Project, and R. V. Welch using methods devised 
by Kuno Foelsch^’ and R. Sauer for the design of axially symmetrical 
supersonic nozzles. 

Following the method of Foelsch, the complete nozzle contour from 
throat to exit was designed for uniform parallel flow at an exit diameter of 
0.5015** and an exit Mach number of 2.8 with isentropic flow. It was esti- 
mated that the boundary layer effect in this nozzle would reduce the 
isentropic exit Mach number to an actual value of 2.6. After manufacture 
and installation of this nozzle in the improved recovery factory apparatus 
the Mach nximbers measured at Station #1, which is 0.81 inches from the 
nozzle exit plane, range from 2.68 at high to 2.33 at low rates of flow 
indicating that the average actual exit Mach number is slightly less than the 



desired value of 2.6. 




t 



The Foelsch method linearizes the characteristic equations of axially 
symmetrical flow by comparing the conditions of flow in a cone with those in 
a nozzle. Mathematical expressions for the nozzle contour can thus be 
obtained. The design begins with an assumed source type flow OA (See 
Fig. E-1) which is altered in region AE to give parallel exit flow. The 
straight conical section AD assures conical source flow on entering 
region AE. The throat is designed by using the equation of continuity to 
determine the throat diameter and constructing TD as a circle tangent to OA 
at point D and parallel to the axis at throat T. Circular arc TC is a con** 
tinuation of circular arc TD, and circular arc Cl is drawn to produce an 
inlet nozzle diameter of 1.250 inches. 

The Sauer method is a part graphical, part analytical, solution of the 
characteristic equations of axially symmetrical flow assuming that the 
flow is isentropic, supersonic and irrotational; it uses successive approxi- 
mations which converge rapidly. Using the Foelsch contour from points T 
to A as a starting point, the exit portion of the contour was designed by the 
Sauer method for an isentropic exit Mach number of 2.8, The contour was 
then adjusted to obtain uniform parallel flow at an exit diameter of 0.5015”. 
A more complete explanation of the Sauer method follows in Appendix F. 

Figure E-2 compares the nozzle contours obtained by the Foelsch 
and Sauer methods and also shows the Klingensmith nozzle used in the 
original recovery factor apparatus. The Sauer contour is about 10% 



shorter than the Foelsch, but the latter has smaller cross-sectional areas 



between the throat and exit planes. Since both contours were obtained by 
assuming isentropic flow in the nozzle, it was felt that the slightly longer, 
but smaller area, Foelsch contour would have less boundary layer effect 
than the shorter, but larger area, Sauer contour. For this reason. 
Professor A. H. Shapiro, Department of Mechanical Engineering, MIT, 
recommended manufacture of the Foelsch nozzle. 

The nozzle was manufactured of stainless steel by the Cummings . 
Machine Works, Boston, Massachusetts, to tolerances of "''.0005 on inside 
diameters and of - .0001 on inside concentricity. 

It is interesting to note that the design of the nozzle contour by the 
Foelsch mathematical computation method required about 15 man-hours 
while that of the Sauer graphical- analytical method required in excess of 



160 man-hours. 




t 



1 



APPENDIX F 



Description of Sauer Method of Nozzle Design 



The purpose of this appendix is to outline the Sauer method described in 
reference 13 which was used in the design of an axially -symmetric, super- 
sonic nozzle with uniform parallel flow at the exit and with an isentropic exit 
Mach number of 2,8. The reference does not fully explain the application of 
this method to this type of nozzle to those unfamiliar with the theory of 
characteristics. Therefore, for the future reference of those who might wish 
to design a similar nozzle, a brief description of the theory of the Sauer 
Method and a more detailed description of its application to this type of 
nozzle is included in this appendix. The nozzle exit contour found by this 
method appears in Appendix E. 

List of symbols ; 

(p - potential function 

X - longitudinal cylindrical coordinate 

r - radial cylindrical coordinate 

c - velocity of sound 

- Mach angle 

- oblique -angular coordinate system tangent to Mach lines at 

each point 

- component of velocity vector parallel to ''p axis = 



W 



- 



4 > 



^ (p 



Mt 



- first approximation of Mach net 
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Mjj - second approximation of Mach’s net 
Gj - first approximation of velocity net 
Gjj - second approximation of velocity net 

- angle which velocity vector makes with axis in velocity net 



Part A -Theory: 



Flow is assumed to he steady, isentropic, supersonic, irrotational, and 
axially -symmetric. 

From the definition of the potential function , equation of continuity, 
the Euler Equation, and the definition of sound velocity, there results the 



potential equation: 

, 2 . 






m 









I- 



Ur; 



z I j_ 

ixdy u r ir Eq. (1) 



Assuming at each point in the physical plane a new coordinate system 
( , n) tangent to the Mach lines at each point, the potential equation is 
transformed to: 



_ sim 

r Eq. (2) 



If in a given region, we consider that this new coordinate system remains 
fixed while a given point P undergoes a small displacement, there results: 

For a finite difference, the above may be written: 



A — 


SI y 

r 




Eq. (4) 


A W^p = 


y 

r 







The significance of these quantities is shown in Figures F-1 and F -Z. 



Consider a Mach net in the physical plane. A constant velocity is assumed 
in each Mach quadrilateral; hence this area will map in a point in the 
velocity plane. The center of each area is assumed to be at the intersection 
of the diagonals. and are the distances between the centers of two 

adjacent areas measured parallel to the 'p) and axis respectively, and are 
considered positive in the direction of flow. V is the component perpendicula 
to the axis of symmetry of the mean velocity vector between the two adjacent 
areas, r is the distance from the axis to the center of the common side of 
two adjacent areas. 

Part B -Application: 

To illustrate the application of this method to the design of a nozzle, 
the steps in the design of a nozzle will be considered in some detail; values 
being taken from the design of this nozzle. 

I. Source flow in a conical nozzle. 

A. First approximation of flow in a conical nozzle: 

From the design by Foelsch method, the conical source flow 
was found to have a half cone angle of 11°26’ and the Mach number 
at point A was 1.87 5, which results in a Mach angle 
= SIh' — = 32 ." ftf' 

For the first approximation in the physical plane, Mj, Figure F-3, the 
Mach lines are drawn parallel to each other with an inclination of i to 
the x-axis. Having this net, and r can be determined. 
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Points 1, 2, and 5 can be immediately located in the velocity plane. 
Figure F-4, since at each of these points the Mach number and the direction 
of flow is known. An initial approximation for point 3 is located by the 
intersection of characteristics through 1 and 2, which are inclined at ^ 
and - respectively to the x-axis. From point 3 thus located, V for 2/3 
and l/3 is determined, and A p for 2/3 and ^ for l/3 are calculated 
from Eq. 4, and point 3 is relocated. 

The initial approximations of subsequent points are determined by 
extending the characteristics of previously determined points and are then 
corrected in the manner explained for point 3. The points so located deter- 
mine the velocity plane Gq. By measuring V and from Gq by means of a 
subsidiary ellipse, the points may be relocated once again to produce 
velocity plane Gj. 

B. Second approximation of flow in a conical nozzle: 

From Gj, the direction of the velocity vectors, 6, at each point 
is measured, and the Mach angle'/ is determined by using a sub- 
sidiary ellipse. The direction of the characteristics at each area 
in the physical plane is given by 6 i'/ . The new Mach net Mjj is 
constructed by drawing the new Mach lines with directions as 
determined by interpolating between the corrected directions. 

After locating the new centers of the quadrilaterals in the Mu 
plane, they are located in the Gj plane to produce the Gj’ net. 

The points of the Gj’ net have the same relative position with 



respect to the points of the Gj net, as do the centers of areas in 
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the Mjj net have with respect to the center^areas in the Mj net. 

To construct Gjj, A W^p and are calculated from Mjj and 

Gi’. The corrected Mach directions found from Gjj are drawn 
through the centers of areas of Mjj. 

C. Third and higher approximations of flow in conical nozzle: 

P'or the third and higher approximations, the procedure for the 
second approximation is continued until for a given approximation 
there is negligible correction necessary to the Mach net in the 
physical plane. 

II. To correct the conical nozzle for parallel outflow: 

The corrected nozzle results from joining a new Mach net M* to the 
Mach net Mjjj (or higher), so as to obtain parallel flow at exit with an 
isentropic Mach number Mj^ = 2.8. The point of isentropic Mach number 
Mj, = 2.8 is located on the x-axis of the velocity plane of Gjjj and from 
corresponding characteristic lines is located in the physical plane of 

Mill, 

To construct Mj*, draw characteristic a-a thru Mj^. (See Figure F-5.) 
The centers of areas at junction of Mjjj and Mj nets are to lie on line a-a. 

In Ml* the directions of the Mach lines are foxmd by extending the character- 
istics of Mjjj. Upstream from a-a, the Mach net and velocity nets remain 
unchanged. Downstream from a-d, the Mach and velocity nets are deter- 
mined by new boundary .conditions; namely, uniform parallel flow of 
Mj^ = 2,8 along a-b. 



A rough estimate of Gq* is made remembering the new boundary- 
conditions. 

Gj* is constructed with values of AW^ p and ^W^p calculated by means 
of Mj* and Gq*. Attention must be paid to the algebraic signs of and . 
Construction must start at point 

With corrected directions of Mach lines in Mj*, the second approxima- 
tions and higher are obtained in a manner similar to that outlined for conical 
flow. 

When the correct Mach net is obtained, the contour of the nozzle is 
drawn. In each area through which the contour passes, the direction is 
determined by the velocity vector for that area. The ratio of exit area to 
throat area (as determined by the Foelsch Method) should be that required 

r 

to give the desired exit Mach number. 
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APPENDIX G 
Nozzle Calibration 

Procedure 

Prior to calibrating the nozzle it was installed in the recovery factor 
apparatus in order that the calibration might yield flow coefficient data which 
included the influence of all components of the apparatus on the mass flow 
rate. Asa result the process was not strictly a “nozzle calibration.” 

In order to determine the flow coefficient, i.e., the ratio of actual to 
isentropic mass rate of flow, twenty-three experimental runs were made 
using both the gasometer^ ^ in the MIT Mechanical Engineering Laboratory 
and a sharp-edged orifice to measure the actual flow rate. This appendix 
will deal only with the data obtained from the gasometer; the results obtained 
from the orifice, described in Appendix J, were not sufficiently consistent to 
be of any immediate use. 

Ten of the rxms, NC-1 through NC-10, were made with apparatus set up 
as indicated in Figure G-1. In these runs high pressure air was hehted in the 
constant temperature bath and passed through the nozzle, test section, and 
orifice to the gasometer, the rate of rise of which determined the actual mass 
flow rate. In the remaining eleven rxms, NC-11 through NC-21, the apparatus 
was set up as indicated in Figure G-2; under these conditions air was passed 
from the gasometer which had previously been filled from an auxiliary air 
line, through the constant temperature bath, nozzle, test section, and orifice 
to a steam-driven ejector. In this case the actual mass flow rate was 
calculated from the rate of fall of the gasometer. 



The latter method more closely simulated the conditions under which 
the recovery factor runs were made since the ejector maintained low 
pressures in the section downstream of the nozzle. For this same reason 
the pressure ratio across the nozzle was always belovy critical and the con- 
ditions for maximum isentropic flow rate always existed. However, since 
the pressure of the air in the gasometer is limited by the number of weights 
which can safely be placed on the tank, the range of nozzle throat Reynolds 
numbers possible was very small and included only values which were 
generally below 10^. The" range of Reynolds numbers was extended some- 
what by using ice in the constant -temperature bath to reduce the temperature 
and viscosity of the air. 

In the former method the lack of a pressure tap in the nozzle makes it 
impossible to determine the isentropic mass rate unless the conditions for 
maximxim isentropic rate exist. As a result, it was necessary to experi- 
mentally determine the upstream stagnation pressure below which supersonic 
flow did not exist at the first station in order to determine the range of 
pressures for which the isentropic flow rate could be computed; it was found 
that the value was 61 psia. Since it is possible to simultaneously have 
subsonic flow at the first station (the first point in the test section after the 
nozzle throat at which there is a pressure tap) and critical flow through the 
nozzle, there is a small range of upstream stagnation pressures for which 
the flow, though maximum isentropic in value, could not be computed. 
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Calculations 

The calculation of the actual mass rate of flow is based on the assump- 
tion of a perfect gas. From the temperature and pressure of the air in the 
gasometer the density can be calculated and from the rate of rise or fall 
the volume rate of flow can be determined since the cross-sectional area of 
the gasometer is fixed at 19.79 square feet. The product of the volume rate 
and density yields the mass rate of flow. 

The maximum isentropic mass rate of flow and the nozzle throat 
Reynolds number are computed by the method outlined in Appendix B. 

A summary of the results of the calculations is tabulated in Table G-1 
and the plot of flow coefficient as a ftmction of nozzle throat Reynolds ntimber 
is shown in Figure G-3. 

Discussion of Results 

From Figure G-3 it can be seen that no calibration data was obtained 
for the Reynolds number range from 1.1 x 10^ to 4.0 x 10^. Although it was 
possible to obtain flow at Reynolds ntimbers in this range, the physical 
limitations of the apparatus made it impossible to calculate the isentropic 
mass rate of flow. As a result, it was necessary to take advantage of the 
data obtained by Ketchum^® in calibrating a similar nozzle. 

In carrying out the calculation of recovery factors the relationship 
between flow coefficient and nozzle throat Reynolds number determined by 
Ketchum was used. The straight line shown in Figure G-3 represents this 



function. 



T 
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The results of the calibration carried out for this report are indicated 
by the points plotted in Figure G-3. It can be seen that the naaximvim deviation 
from the Ketchum line of any point is roughly 1.1% and that the mean deviation is 
considerably less. For this reason and in view of the fact that no better data 
was available at the time of the calculations, it is felt that the use of Ketchum*s 
data is justifiable. 

Recommendation 

In order to determine the flow coefficient as a function of nozzle throat 
Reynolds number over the range for which no data was obtainable by means 
of the gasometer, the following procedures are possible: 

1. A pressure tap can be located in the nozzle such that the isentropic 
mass rate of flowjmay be determined when the flow is not critical; or 

2. The nozzle can be removed from the recovery factor apparatus and 
calibrated independently of the apparatus. 

In view of the fact that it is undesirable to disturb the inner surface of 
the nozzle, it is recommended that the second procedure be carried out to 



calibrate the nozzle. 
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Summary of Results of Nozzle Calibration 
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APPENDIX H 



Thermocouple Calibration 

Thermocouples used for the measurement of the upstream and down- 
stream stagnation temperatures and were those previously 

2 

calibrated by Junge and Margolskee . Results of this calibration are also 
found in Appendix C. 

Thermocouples for use in the lucite test section of the improved 
recovery factor apparatus to measure station temperatures and station 
radial temperature gradients were made of the same spools of #30 con- 
stantan and #30 copper wire as were used for the T and T , thermo- 
couples. The manner in which the lucite test section thermocouples were 
made and installed in the apparatus follows that reconrunended by 

4 

Shovilberg . The station temperature thermocouples were calibrated in 

the MIT Heat Measurement Laboratory with the hot junction in a steam 

bath and the cold junction in an ice bath. The method of calibration follows; 

Symbols ; p^ - barometric pressure (corrected) 

Tp - temperature of boiling point of water °C 

mv - average potentiometer millivolt readings of the 
eleven stations while at steam temperature 

Tj.p - Bureau of Standards Conversion Table (for #30 
constantan and #30 copper wire thermocouples) 
temperature corresponding to w'Vj. , i.e., 

average temperature indicated by the eleven 
stations while at steam temperature 

mv^ - potentiometer millivolt reading at any tempera- 
ture between 0 and 100°C. 



tj - difference between and 
d rp p 

tj. - temperature correction to be applied to t^. 

tj. - Bureau of Standards Conversion Table 
temperature corresponding to mv^ 

t^ - true temperature at the station = tj. --t^ 
p^ = 758.292 mm Hg 

Tp = 100,0 +0.0369 {P|j-760) -0.000020 (pb’760)2 = 99.937 °C 

mVj,p = 4.2884 mv 

Typ = 100.27 °C 

t , = T - T = 0.33 °C 
d rp p 

Assuming (a) that at 0 °C, the potentiometer reading (mv^) will be 
zero so that = 0 and 

(b) that the temperature correction t is a linear fvmction 
of mVj,, then 

V td 0.33 °C °F 

^ = =___ = .1385 

mv_ mv«_ 4.2884 mv mv 

T rp 

therefore ^ ~ ~ ^c ~ ^r ~ naVj, (.1385) Op 

A table of values of t^ vs mv^ was then computed by assuming 
values of mv^ at 0.1 mv intervals, obtaining t^ from the Bureau of Standards 
Conversion Table and then using the above equation to compute t^. (See 
Appendix C.) Linear interpolation between these values results in errors 
much smaller than the accuracy of the thermocouple measurements. 

At steam temperature the average deviation between the station 



millivolt reading and the averaged millivolt readings of the eleven 



I stations (mv^p) is i 0.0025 mv or i 0.09 ^F. At temperatures between 
0 and lOO^C, the average deviation will be proportionately less. 
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APPENDIX I 

Vacuum System Precautions 

The vacuum system itself is adequately described in Chapter V. The 
purpose of this appendix is to give the more important precautions to be 
observed in the operation of the vacuum system in order to maintain a low 
vacvium of 0,003 to 0.005 mm Hg pressure. 

One of the major sources of potential trouble is the diffusion pump 
primarily because the pump is made of metal and its operating fluid cannot 
be observed nor can mal-operation of the pump quickly be detected by 
means of the vacuum measuring system. The condition and amount of 
octoil in the pump is not readily determinable since inspection requires 
breaking the seal on the not too accessible pump drain plug, removing the 
octoil and replacing the oil after measuring and inspecting. Inspection of 
the trap placed between the diffusion pump and the megavac pump to 
prevent carry over of octoil into the megavac pump is also difficult since 
its drain plug is quite inaccessible. Because of the danger of disturbing 
the joint seals in nearby piping and also because about 5 hours of running 
time is required to remove all dissolved gases from any fresh octoil placed 
in the system* removal of the octoil from the diffusion pump for inspection 
is not recommended at intervals of less than two months unless some of 
the following precautions have been violated or the pump is definitely not 
operating satisfactorily. Because of the fractionating principle, however, 
one filling of octoil should last indefinitely if not subjected to abuse or 
action of chemical vapors. Therefore, the diffusion pump should be quite 










mh 



f 
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dependable if the following precautions for the prevention of decomposition 
and contamination of the octoil are observed; 

1. Do not turn on the diffusion pvimp heater until the mega vac pump 

t 

has reduced the vacuum system pressure to less than 0,1 mm Hg 
because the hot oil vapors will decompose on exposure to air at 
pressures above 0.1 mm Hg. 

2. In order to condense the hot oil vapors before coming into contact 
with the air, continue circulating cooling water after turning off the 
pxrnip heater for at least 15 minutes before opening the vacuum system 
to the atmosphere. 

3. Open vacuum system to the atmosphere very slowly. 

4. When system is open to the atmosphere, the circulating water must 
not be turned on since water vapor will condense inside the pump and 
contaminate the oil. 

5. It is preferable to leave the vacuum system at pressures below 
atmospheric when system is not in use than to leave system open to 
the atmosphere. 

6. Regulate the circulating water such that the cooling coil and entire 
surface of the diffusion pump beneath the coil are quite cold during 
operation. Water leaving the cooling coil should be slightly warm to 
the touch. 

7. When octoil is to be removed for storage during long periods of 
pump idleness, it should be stored in glass stoppered, amber 



colored bottles. 
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8. Cleaning of the punnp should be undertaken whenever the vacuum 
system is dismantled. Acetone should be used for this purpose, 
followed by air circulation, and then evacuation of the pump. Carbon 
tetrachloride must not be used since this solvent clings to the metal 
parts and is difficult to remove. 

The mega vac pump requires no special care other than proper lubrica- 
tion of the motor, maintaining correct level of clean high vacuum oil in the 
pump, prevention of overheating, etc. The pump should give trouble-free 
and efficient operation for about 100,000 hours. 

The remaining causes of poor vacuum come under the heading of leaks 
through the numerous joints and seals in the vacuum system. Constant 
vigilance and application of glyptal and air -drying varnish to suspected 
joints is mandatory. It is preferable to apply the glyptal and varnish when 
the vacuum system is open to the atmosphere and time should be allowed for 
thorough drying before placing the system in operation, if possible. The 
covering of glyptal and varnish should be removed by means of glyptal 
thinner or benzene before dismantling a joint. When assembling any part 
of the vacuum system piping or test section, it is essential to use every 
means possible, such as, cleaning of parts before assembly, use of high- 
vacuum grease on all rubber gaskets and screw threads, proper tightening 
of bolts and compression fittings, etc., to prevent future leakage. 

Some precautions to be observed in the vacuum measuring system are: 
1. Even though a safety trap is provided, care should be taken when 



comparing the pressure in the absolute pressure manometer manifold 
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with that in the vacuixm system to prevent silicone oil being forced into 
the recovery factor test section as a result of the pressure in the manom- 
eter manifold being much greater than that in the vacuum system, 

2. After measuring the pressure in the vacuum system by means of 
the McLeod gage, lower the mercury reservoir to its lower position to 
facilitate evacuation of that pai^t of the apparatus. 

3. The mercury reservoir mentioned above should be lowered slowly 
to prevent an air bubble entering the mercury column from the mercury 
column air cushion and thus raise the pressure in the vacuum system 
appreciably. 



APPENDIX J 



Orifice 



Design and installation 

A standard orifice plate for measuring the downstream flow rate was 
designed and installed in the downstream section of the recovery factor 
apparatus in order to provide a comparison with the calculated nozzle flow 
rate and to detect leaks in the apparatus* 

The orifice was designed in accordance with reference 12 for maximum 
possible flow rate occurring at the highest expected upstream stagnation 
pressure and the lowest expected upstream stagnation temperature. Thus: 



pressure taps was determined for the above flow rate from the following 
equation taken from Section VII of reference 12: 



in which the orifice upstream pressure (P^ q) and the pressure differential 



eters to be used for their measurement. From the area (A^) thus determined, 
the diameter (D 2 ) was found to be about 1.5**. 



Wo = .5319 A* ^us = .5319 (.05662) = .1355 Ib/sec 




The area of the orifice to be installed in the 2 l/2*’-pipe with radius 




( 1 ) 



(APq) across the orifice are selected from a consideration of the manom- 



The orifice plate was then machined from l/8” stainless steel plate 
with an orifice diameter of 1.498 inches and with a 45°-bevel on the down- 
stream face. The plate was installed between two 2 l/2-inch pipe flanges 






with a rubber gasket on each side of the plate. The length of straight piping 
provided upstream of the orifice was 22 pipe diameters while that down- 
stream was 10 pipe diameters, both of which are greater than the minimum 
required by Section IV of reference 12. The. upstream pressure tap was 
installed at a distance of 1 pipe diameter from the orifice plate and the 
downstream tap at a distance of l/2 pipe diameter. 



Calibration 



The orifice was calibrated after installation, using air flow measured 
by the gasometer at the same time the nozzle was calibrated. The orifice 
flow coefficient was calculated from the following equation: 

W^asomeber 



Ko = 



0.668 F), E Y I f?. aR 



(2) 



R Xos 



Results of the calibration are shown in Figure J-1, but it should be noted 



that this calibration is valid only for 



AP, 



PfO 



0.2 which corresponds roughly 



to upstream stagnation pressures (P^g) greater than 25 psia when T^g is 
near 110 ®F. 



Calculation of Flow 

The method of calculating the orifice flow rate is adequately described 
on page 56, reference 12. There remains only the necessity of indicating 
how the values were measured. Referring to equation (1): Pj^q was obtained 
from the orifice upstream pressure tap reading converted to psia; A Pq 
obtained directly from the pressure differential manometer across the 

Cj 

orifice; T^ was assumed to be equal to T(jgJ E, the are multiplier for 



thermal expansion of the primary element due to the pipe temperature, is 
equal to 1 at room temperature; and Y, the expansion factor, was calculated, 
knowing the ratio of orifice to pipe diameter and the ratio of ^ Pq to 
Results: 

As can be seen from Figure J-1, the orifice flow coefficients determined 
from the gasometer flow vary as much as i 5% from the mean value of the 
orifice coefficient. Assuming the gasometer flow to be correct, because it 
yielded nozzle flow coefficients that were quite close to values previously 
determined, the orifice flow determination is not accurate enough to permit 

the orifice to be used for the purpose for which it was designed. 

A P 

The orifice is limited to flow rates in which ^ < 0.2 and thus for 

Pi " 

Pus gi’eater than 25 psia when 110 This limitation can be improved 

by redesigning the orifice for the smaller flow rates. However, it was felt 
that re-design with a smaller orifice area would not improve the large 
variation of the orifice flow coefficient to within the i 0.5% necessary to 
make the orifice usable as a means of determining leaks or errors in the 
determination of nozzle flow rates. Therefore, re-design was not attempted. 
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APPENDIX K 



Modification to the Recovery Factor Apparatus 

After completion of runs F-1 to 35, the recovery factor test section was 
modified to reduce the heat conduction and radiation between the l/S-inch O.D. 

r 

copper tubing used for measuring the station wall pressures, and the lucite 
test section. One end of the copper tubing is exposed to room temperature, 
while the other is fastened to the lucite test section. As discussed in 
Chapter VI, experinoents with dry ice placed on the exposed ends of the copper 
tubing showed that an appreciable amount of heat was being transferred from 
the room to the test section via the copper tubing; or vice versa, depending 
upon the relation between the station wall temperature and room temperature. 

The modification consisted primarily of removing all copper tubing 
inside the recovery factor apparatus, leaving only a length of 1 l/4 inches at 
the upstream end plate and at the station boss compression fitting. The 
removed copper tubing was replaced with low conductivity, l/8-inch plastic 
tubing. The ends of the plastic tubing were slipped over the short lengths 
of the remaining copper tubing and a tight joint was effected by means of 
small tube clamps and a coating of glyptal. This change effectively prevents 
all but small amounts of heat transfer by conduction and radiation between 
the room and the test section by the pressure tubing. 

To reduce the other possible radiation effects in the apparatus, such as 
that between the lucite test section and the textolite outer tube, aluminum 
foil was wrapped loosely around the test section, around each pressure tube, 
and then around the assembly of the pressure tubing and test section. Both 
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end plates were covered with several layers of aluminum foil and a layer 
was loosely secured to the inside of the textolite outer tube. The recovery 
factor apparatus was wrapped with aluminum foil externally, as well. 

In order to determine the temperature surroxinding the lucite test 
section, and also to determine the temperature gradient existing between the 
test section and the textolite outer tube, two thermocouples were installed 
inside the apparatus. The thermocouple wires were led out of the apparatus 
through four 45° holes drilled in the downstream end plate in a manner 
similar to the thermocouple holes drilled in the upstream end plate, and 
thence through the joint between the downstream end plate and the textolite 
outer tube. A tight seal was effected by means of high vacuum grease 
applied to the wires, the end plate gasket and the joint between the end plate 
and outer tube. Externally, this joint was coated with glyptal eind air -drying 
varnish. 

One of the above thermocouples, called Tj^, was located about l/Z inch 
from the lucite test section between the bosses of stations #5 and 6. Two 
layers of aluminum foil separate the thermocouple from the lucite itself. 

The other thermocouple, T-p, was placed on the same horizontal plane and 
in the same longitudinal position as Tj_^ but was fastened by adhesive tape 
directly to the inside of the textolite outer tube. 

After the modification was completed, run F -36 was made at Tus = 
109.89°F and the recovery factors obtained were found to agree quite well 
with those previously obtained at similar upstream stagnation temperatures. 
This indicates that the modification does not appreciably change the data on 



recovery factors already obtained at T^g near 110 ®F. This was the expected 
result since at this T^g the wall temperature in the lucite test section is 
approximately equal to room temperature and thus only small heat transfer 
was possible via the copper tubing. At other stagnation temperatures, however 
the heat transfer before the copper tubing was removed was proportionately 
greater and this modification will result in much more accurate recovery 



factors. 
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APPENDIX L 

Effects of Measurement Errors 

The equipment used to measure observed quantities in the recovery 
factor and friction factor runs is the same as that used by previous investi- 
gators working on other phases of this project. Since the limitations on the 
precision of measurements have been discussed in references 2 and 4, it is 
the purpose of this appendix to present an analysis of the effects of the 
measurement errors on the calculated quantities. 

The recovery factor is dependent on the Mach number, which is a 
computed quantity, and the adiabatic wall temperature and upstream stagna- 
tion temperature, both of which are measured quantities. The Mach number 
in turn is a fxinction of the local wall pressure, upstream stagnation pressure, 
and flow coefficient. 

Consider the equation (See Appendix B): 

Using the technique of taking the logarithm of both sides of the equation and 
differentiating leads to; 

d CPw/Po>s) dCw _ dM t (K- OM^ dM 

Pw/p^i Cw M 2 ^ M 

for a/a* and k both constant. 

For small changes in variables the differentials in the above expression 
may be replaced by delta -qviantities to indicate small finite differences. To 






study the effects of errors in flow coefficient only assume that there are no 
errors in the pressure measurements. Then 

K- \ \ ^ M 

C w ~ ^ Si t- fill W\ 

Assuming k = 1.400 and rearranging; 

^ M _ ( > -V- . Z W\^) A C yy 

(» -k- .4 M^) Cw 

The above equation expresses the interdependence of the percent error in 
calculated Mach number, percent error in flow coefficient, and Mach number; 
the function is plotted in Figure L-1. 

From the same original equation and in a similar manner the effects of 
errors in measurements of wall pressure and upstream stagnation pressure 
can be determined. The resulting error equations are: 

A M = (\ + A 

M ( \ +• . 4 M P^s 

A M ^ ^ ,^M^) A Pyy 

(\ + .4M^) Pw 

These functions are also plotted in Figure L-1. 

To determine the effect of errors in Mach number on calculated recovery 
factors, the same technique is applied to the definition of recovery factor 

< - CTo. /To.s') 



R' 



Since is a function of Mach number and since T^^/T^g 

appreciably constant, the following can be derived: 



IS 




\ - C ' - Tav./To.sH ^ oj 



For T-^^/T^g - 0,94 

A V-<' ^ ( 0-^S9 N 

Figure L.-2 shows the percent error in recovery factor as a function of 
percent error in Mach number for the above value of T^^/T^g* 

From the definition of recovery factor it can be shown that 

^ aT^W 

R' V - I 

/T \jki^ 

This fvmction is plotted in Figure L-3 for various values of Mach number 
and for Taw/Tus constant at 0.94. 

As a result of this analysis the probable errors in calculated recovery 



factors due to the limitations on the precision of measurements are tabulated 
below; 



Measured 

Quantity 


Probable 

Measurement 

Error 


Maximum 
Probable 
Error (%) 


Maximum 
Probable 
R’ Error 


Pus (Bourdon gage) 


0.5 psi 


1.3% 


M = 1 
1% 


M = 2.: 
0.1% 


P (He manometers) 

us 


0. 1 0 cm 


0.2% 


0.2% 


0.05% 


(Hg manometers) 


0, 1 0 cm 


2% 


1.6% 


0.3% 


P^ (absolute manometers) 


0, 1 0 cm 


0,5% 


0.4% 


0.05% 




0.01 


1% 


0.8% 


0.1% 


^aw ('^aw'/'^us “ *94) 


0.4°F 


0.1% 


0,8% 


0.2% 
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